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For a system of N spins 1/2 there are quantum states that can encode a direction in an intrinsic 
way. Information on this direction can later be decoded by means of a quantum measurement. We 
present here the optimal encoding and decoding procedure using the fidelity as a figure of merit. We 
compute the maximal fidelity and prove that it is directly related to the largest zeroes of the Legendre 
and Jacobi polynomials. We show that this maximal fidelity approaches unity quadratically in 1/N. 
We also discuss this result in terms of the dimension of the encoding Hilbert space. 
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I. INTRODUCTION 



Entanglement and superposition are the most characteristic features of quantum states. They play a central role 
in the storage and transmission of information in the quantum world and are responsible of the many remarkable, 
and often intriguing, quantum effects that are constantly being discovered. These effects, in turn, are providing new 
insights in the difficult task of understanding quantum information. 

Some time ago Peres an Wooters Q posed an interesting question. Imagine a quantum system composed of 
several subsystems, which are not necessarily entangled. How can we learn more about this system? by performing 
measurements on the individual subsystems or on the system as a whole? They showed evidence that the latter, the 
so-called collective measurements, are more informative. Obviously entanglement is the property responsible for this. 
In this case, however, it is not explicit, since the system can be chosen to be in a product state, but hidden in the 
collective measurement. 

Later Massar and Popescu || addressed a more concrete problem. Imagine Alice has a system of N parallel spins. 
She can use this system to tell Bob the direction along which some given unit vector n is pointing. She just have to 
rotate, or prepare in some other way, the state of her system so that it becomes an eigenstate of ft ■ S, the projection of 
the total spin in the n direction. The state is then sent to Bob whose task is to determine the direction encoded in the 
state. He will need to perform a collective measurement and from each one of its outcomes, labeled with an index r, 
he will have a guess for Alice's direction given by a unit vector n r . To quantify the quality of Bob's measurement 
Massar and Popescu used the average fidelity which is defined by F = J2r I dn(l + n ■ n r )/2 P r (n), where Alice's 
directions n are assumed to come from an isotropic source. Here P r (n) is the probability of getting the outcome r if 
Alice's direction is ft, and dn is the rotationally invariant measure on the unit 2-sphere. The authors proved that the 
maximal average fidelity Bob can achieve is F — (N + 1)/(N + 2) which is readily seen to approach unity linearly: 
F ~ 1 — 1/N. Explicit realizations of the optimal measurements with a finite number of outcomes were obtained in (3] 
for arbitrary N and minimal versions of these measurements for N up to seven are in Q] . 

A new surprise was recently presented in In this paper the authors consider N — 2 and show that states with 
two anti-parallel spins, | fj.), | J,f) provide a better encoding of Alice's directions than the two parallel spin states used 
in ||, ||, ||. The average fidelity is now (3 + y/3) /6 which is larger than 3/4 for two parallel spins, i.e., Bob can have a 
better determination of Alice's direction if she uses anti-parallel spins. This is a startling result, since classically one 
would expect that a direction is encoded equally as well in a state pointing one way as in one pointing the opposite 
way. The main reason why this is not so in the quantum world, as it will become clear from our work, is the different 
dimensionality of the Hilbert spaces to which two parallel or two anti-parallel spin states belong . 

At this point, the obvious reaction is to ask ourselves what are the best states Alice can use to encode directions. 
Since the very natural state with only parallel spins is not optimal for N = 2, we expect that neither will it be for 
arbitrary N. Hence, one has to search the optimal encoding state among all the eigenstates of ft ■ S that Alice can 
produce. In this paper we present this very general analysis. We compute the maximal average fidelity (hereafter we 
will usually drop 'average' when there is no ambiguity) for arbitrary N and show that it approaches unity quadratically 
in 1 /N, as compared to the linear behavior found in Q for parallel spins. As a byproduct, we also compute the maximal 
fidelity for encoding states of two arbitrary spins s such as two nuclei. A short description of the main results of this 
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analysis was presented in || . These results have been recently corroborated by numerical analysis . 

The paper is organized as follows. In section II we introduce our notation and conventions and present a detailed 
calculation of the maximal fidelity for TV = 2. We show that the fidelity obtained by Gisin and Popescu in (s| is 
optimal (a result also obtained in || using different methods) . In section III we analyze the more general case of two 
spin s states. The analysis for any number N of spins is in section ^ and our results and discussion are in section 
We conclude with an appendix containing technical details. 



II. TWO SPINS 



We start by assuming that Alice has two spins in a general eigenstate of n ■ S (We skip the analysis of the simplest 
situation in which Alice has only a spin. The reader can find it in ^| , and our general formulae of section [iy| can 
be also particularized to this case). We can think of it as a fixed eigenstate of S z = z '• S (z is the unit vector pointing 
along the z direction) that Alice has rotated into the direction n — (cos 0sin0, sin ^>sin0, cos 8). It is convenient to 
work in the irreducible representations of SU(2). In the present case, l/2®l/2 = 1®0, the general form of this 
fixed eigenstate is 

\A) = A+]l,l> + A)|1,0) +A-\1,-1) +A a \0,0), (1) 

where, as usual, the normalized states of the basis, \j, m), are labeled by the total spin S and the third component 
S z : S 2 \j, m) — j(j + to) and S z \j, to) = m\j, to). In the following we stick to the general form ([!]) to treat all the 
cases jointly, but one should keep in mind that only combinations with definite S z will be relevant for our analysis. 
The rotated state U(n)\A), where U(n) is the element of the SU{2) group associated to the rotation z — > n = R z, is 
precisely Alice's general eigenstate of n-S. Obviously U(n) is reducible since it has the form U(n) — (n)(BU^ (n), 
where denotes the SU(2) irreducible representation of spin j. 

Next, Alice sends the rotated state to Bob, who tries to determine ft from his measurements. The most general 
one he can perform is a positive operator valued measurement (POVM). We specify this POVM by giving a set of 
positive Hermitian operators {O r }, that are a resolution of the identity 

I = Y,Or. (2) 

r 

For each outcome, r, Bob makes a guess, ft r , for the direction. As we brought up in the introduction, the quality of 
the guess is quantified in terms of the fidelity which we can view as a 'score'. To Bob's guess n r , we give the 'score' 
/ = (1 + n ■ ri r ) /2. We see that he fidelity / is unity if Bob's guess coincides with Alice's direction and it is zero when 
they are opposite. Thus, if n is isotropically distributed the average fidelity can be written as 



, 1 + ft ■ ftr , , ^, „ , 

dn r -tr[p(n)O r ], (3) 



where pin) = U(n)\A)(A\U' (n) and dn was defined in the introduction. The evaluation of F can be greatly simplified 
by exploiting the rotational invariance of the integral (0). If we define R r through the relation 



and make the change of variables 



We have 



ri r = R r z (4) 

R^n^ft. (5) 
1 + ft ■ z 



r J 

where 

tt r = U\n r )O r U{fi r ). (7) 

Notice that in general ^ r Vt r ^ I. We can regard f2 r as fixed or reference projectors associated to the single 
direction z. In this sense, they are the counterpart of Alice's fixed state \A). Inserting four times the closure relation 
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^ k \k)(k\ = I, where k — +, 0, — , s, and {\k}} is the basis of the representations 1 ® 0, 



we obtain 



l±> = |1,±1) 
|0) = |1,0) 
\s) = |0,0), 



F = g A* At u k3 Jdn 1 + ™ Sd ®* kl {n)^ji (n) 



(8) 



(9) 



Here the indices fc, i, j, I, run also over +,0,— ,s; ®kj(n) = [D^ © D^]kj(n) — (k\U(ri)\j) are the SU(2) rotation 
matrices in the 10 representations, and 



= ^(fe|n r |i). 



(10) 



Now, one can easily evaluate the integrals and obtain the fidelity 

F = A^WA, 

where A = (A + , Aq, A_, A s Y and A''' is its transposed complex conjugate. The matrix W is 



(11) 



W 



/ 3w ++ + 2^oo + w_ 
12 



V 



U)++ + OJQQ + 0J— 

6 



* 

2cjoo + 3w_ 



12 



* 


\ 






IT 




* 
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(12) 



where the entries marked with * are not relevant for our analysis since we only consider eigenstates of S z for the fixed 
states | A) . These, and the corresponding rotated states U(n) \ A) , are the only ones that point along a definite direction 
in an absolute sense, i.e., even if Alice and Bob do not share a common reference frame. ^From its definition jic|), 
it follows that u>jj are real nonnegative numbers but u>ij are in general complex numbers for i ^ j. There are other 
constrains on u>ij stemming from the condition Y) O r = I: 



w ss = 1, ^ = 3. 

i=+,o,- 



(13) 



Because of the Schwarz inequality, we also have 



|^ 0s | 2 < uJon^ss = woo- (14) 
Let us discuss the implications of these equations for different values of m. 

Case m — ±1 

The fixed state \A) for m — 1 is simply \A) = |1, 1), i.e., A + — 1 and Ao = A- = A s = 0. In this case the fidelity 
is given by the element W ++ of (12), 



F = W ++ = 



3w ++ + 2w o + w_ 



3 t^oo + 2cj_ 

4 ~ ' 



3 

^4' 



12 4 12 

where the second condition in ( |l3| ) has been used. The maximal value, which we denote by F+, is then 

+ 4 

This value occurs for 



(15) 



(16) 



u — = wqo = 



oj ++ = 3. 



(17) 
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The case m = — 1, for which \A) = |1, —1), is completely analogous with the index substitution + <-> — . The maximal 
value of the fidelity is also F_ =3/4. 

Case m — 

For m = one has \A) = Ao|l,0) + A S \Q, 0), with \Ao\ 2 + \A S \ 2 = 1. The maximal fidelity is the largest eigenvalue 
of the 2x2 submatrix of ([L2h corresponding to the m = subspace: 



(18) 

It reaches its maximal value, F , for 

ujoq = 3 w++ = = 0. (19) 



Substituting back into (|l^) we obtain 



The corresponding eigenvector is 



(20) 



1 e 



i<5 



|A) = -^|1,0) + -^|0,0), (21) 

where the phase is the unconstrained parameter S = argw s o Notice that the family of states (|2l]) contains entangled 
as well as unentangled states. With the choice e lb = ±1 one obtains the product states | fj.), | J.T); precisely those 
considered by Gisin and Popescu || , which led them to the conclusion that anti-parallel spins are better than parallel 
spins for encoding a direction. 

^From this analysis one can also obtain important information about the optimal POVM. Taking into account that 
one can always take the projectors O r to be one-dimensional ||, we can write Bob's reference projectors f2 r as 

fi r = c r |$V)(* r |, (22) 

where |>t r ) are normalized states and c r are positive numbers. The values of u>ij (see Eq. |l^) endow the information 
about the components of |^ r ) in the spherical basis (||). To be specific, consider states with m = 0. The maximal 
fidelity condition ( [l9| ) implies that the states |\& r ) must have also m = 0, hence |^ r ) = a r |l,0) + /3 r |0,0). This result 
is, to some extent, what one expects: In order for a POVM to be optimal, the measurement must project on states as 
similar as possible to the signal state. Further, the Schwarz inequality (|l4|) becomes equality if and only if a r j (3 r = A 
for all r. If this is the case, the fidelity can reach the maximal value Fq. Then, imposing the POVM conditions dl3| ) 
it is straightforward to verify that all \^ r ) must coincide with a single state, which we denote by \B), 

\9 r ) = \B) = ^\l, 0) + ^|0,0) ) (23) 

The relative weights of the |1, 0) and |0, 0) components, v3 : 1, are easily understood as being the square root of the 
dimension of the Hilbert spaces corresponding to j = 1 and j = 0. We therefore see that optimal POVMs can be 
obtained by rotating the single reference state \B). The weights c r are free parameters except for the constrain 



4. (24) 



Because the Hilbert space has dimension four, a POVM (optimal or not) must consist of at least four projectors. 
Let us show that indeed an optimal POVM with this minimal number of projectors exists. Since the number of 
projectors in the POVM equals the dimension of the Hilbert space, we are actually dealing with a von Neumann 
measurement, i.e., 

O r O s = O r 6 rs . (25) 

Hence, (^1^) = 1 => c r = 1 for the four values of r, which is, of course, consistent with (24). Inverting (Q) and 
taking into account (p2j), we see that the four unit vectors n r have to be chosen so that 

4 4 

^2o r = ^2u{n r )\B){B\U^(H r ) = I. (26) 

r—t r—1 
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By symmetry, they should correspond to the vertices of a tetrahedron inscribed in a unit sphere, i.e., n r — 
(cos 4> r sin 9 r , sin <p r sin r , cos 9 r ) with 

cos 0i = 1, 0i = 0; , , 

cos0 r = -§, <t> r = (r-2)^, r = 2,3,4. 1 u 

It is easy to verify that with this choice condition ( |26| ) is fulfilled and the maximal fidelity (|2(j) is attained. One can 
check that the four projectors ( ^6|) are equal to those already considered by Gisin and Popescu in ||. Our aim here 
was just to present a motivated explanation for their choice of POVM. Finite optimal POVMs for N > 2 are less 
straightforward to obtain. However, the results of j|, which enables us to construct finite POVMs for code states 
with maximal m, \N/2, N/2) = | f f . N .\ f), can also be used here for other values of m. We will comment on this issue 
in our last section. 



After dwelling on minimal POVMs, it is convenient to consider also the other end of the spectrum: POVMs with 
infinitely many outcomes or continuous POVMs juj. They will be used in the general analysis in the sections below, 
where they will prove very efficient. Recall that for any finite measurement on isotropic distributions it is always 
possible to find a continuous POVM that gives the same fidelity ||. Therefore, restricting ourselves to this type of 
measurements do not imply any loss of generality. We illustrate this point for N = 2 and m = to introduce the 
notation that will be used in the following sections. 

We have seen that the matrix elements w, 3 contain all the information required for computing the fidelity, indepen- 
dently of any particular choice of POVM. Any measurement for which uiij satisfy the condition ( |l7| ) for m = 1 or jl9| ) 
for m — is surely optimal. A continuous POVM is just a particularly simple and useful realization. It amounts to 
taking the index r to be continuous, i.e., 

r J 

where the subindex B in the invariant measure refers to Bob (measuring device). Substituting (|2^) into ( |Toj ) one 
obtains in the continuous version 



Ukj = J dn B c(n B )(k\B)(B\j), (29) 

where \B) is the normalized state (|2^) and c(ng) is a continuous positive weight, which plays the role of c r and 
according to ( pi]) must satisfy 

J dn B c{n B ) = 4. (30) 
We now show that in fact c(n B ) is a constant and, hence, equal to 4. Condition ( p6| ) reads 

dn B c(n B )U(n B )\B){B\U^(n B ) = I, (31) 



which is equivalent to 



2j + l 



dn B c(n B ) T) ( £ (n B ) T) ( £l*(n B ) = 5 jf 6 mm ,, j,f =0,1. (32) 



4 

Using the well known orthogonality relation of the matrix representations of SU(2) [tl|, 



dn 2)W m2 (n) S)^; 2 (n) = -L- 5jJ ,S mim , : (33) 



c(n B ) = c = 4, (34) 



one obtains 



which is just the total dimension (3 + 1) of the Hilbert space to which the state (|23j) belongs. Therefore, the projectors 
0(n B ) — cU{n B )\B){B\U\n B ) in ( pi] ) describe an optimal continuous POVM. They are obtained from the fixed 
state (||^) in a manner analogous to the construction of the minimal POVM in (^6|) and (p7|), excepting the constant 
factor c required by the normalization of the matrix representations of SU(2). 
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To complete the analysis of TV = 2, we calculate the maximal fidelity for a given (non-optimal) fixed state \A) 
with m = 0. Without any loss of generality it can be written as 



\A) = \A \ |l,0) + |A s |e* 5 |0,0), \A \ 2 + \A S \ 2 = 1- (35) 

where we have used the same phase convention as in (21). From (|l2|), and the constrains (|l^) and (p^), it is 
straightforward to see that the maximal value of the fidelity is 

f4 = i + ^H (36) 

To attain this value, Bob must perform an optimal POVM, characterized by (p3|). He may use, for instance, the 
minimal one (Eqs. |26|-p7|), or the continuous one, 0(tib)- This result shows that for any fixed state ( |35| ) with 
1/2 < \Aq\ < v / 3/2 the fidelity is higher than that of the parallel case (i.e., m = ±1) for which F = F± = 3/4. 

III. TWO SPINS S 

Imagine now that instead of two spin- 1/2 Alice can use two equal arbitrary spins si = S2 = s to encode the 
directions. This can be seen as a generalization of the simple case studied in the preceding section. However the most 
important feature of this analysis, as it will be shown in section [W, is that it provides the solution of our original 
problem, namely, that of obtaining the maximal fidelity when Alice has N spin- 1/2 particles at her disposal. 

According to the Clebsch-Gordan decomposition, a normalized eigenvector of the total spin in the z-direction with 
eigenvalue can be written as 



\a)= Mwa); £ l^l 2 = 1 ; ( 37 ) 



j=rriA 3=mA 



where J = 2s. The state \A) and its components Aj should carry the label ttia to denote the different eigenvalues of 
S z , however, we will drop it to simplify the notation. A general eigenstate of ft ■ S has the form U(ft)\A), where U(n) 
is now 

U{n)= U^(ft) (38) 

j=m A 

The POVM projectors can be constructed from a fixed state \B) of the form 

.7 

\B)= £ Bj^ms), (39) 

j=m B 

namely, O(ns) = cU{fiB)\B){B\W (Hb)- Note that \B) is an eigenvector of S z with eigenvalue ms, though we also 
drop the label ms here. The absolute value of the coefficients Bj and the positive weight c are determined by the 
completeness relation f driB 0{ub) = I, which using ( |33"| ) leads to the normalization condition 



|2*il = y^, (40) 

and a value for c given by 

c=(J+l) 2 -m|. (41) 



Notice that the factor 2j + 1 in ( |40[ ) is just the dimension of the Hilbert space of the irreducible representation j of 
SU(2), and c is the dimension of the total Hilbert space. Thus, ( |39| ) is the straight generalization of the states (pS))- 
The fidelity can be written as 



F = c £ 

jj'=m 



A,A),B]B y J d nl±^^ Bm Jn)^* mA (n), (42) 
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where 



m = max(m J 4, mj). 



(43) 



The integral in ( |42| ) can be easily computed by noticing that cos 6 — D^^n). Using again the orthogonality rela- 
tions we have 



dn cos 



0D« m2 (n) D^ 2 (n) = ^7<10; jmi|i'mi><10; jm 2 |j'm 2 >, 



were (jirai; 1 J3"t-3) are the Clebsch-Gordan coefficients of j± (g> ji 2 — > jz- The fidelity can be cast as 



1 1 



F = 2 + 2 £ Mj|Aj|2 + ^ ^ ( A *j-l A rf + A i-i A >i 

j=m j= m+1 

where the last term is zero for uia < wis and the coefficients fij and i/j are 

mAVflB 



1 jn— 1 



J=m A 



2 



flj 



j(j + 1) 

e l5j ( ( j 2 - m 2 A )(j 2 - m%) 



1/2 



(44) 

(45) 

(46) 
(47) 



The phases Sj in (|47|) are arbitrary. They are just the generalization of the single free phase of (f23J). Here we have 



arg(£>*i?j_i). The maximal fidelity is achieved by choosing 8j equal to the phases of the signal state \A): 



5, = argOB*^) = arg^A^) 



(48) 



We see now that all terms in (^5|) are explicitly positive with the exception of the last one, which necessarily vanishes 
for optimal states \A), i.e., Aj = for j < m. Gathering all these results, we obtain for the fidelity: 



F = - + - A* MA. 
2 2 

Here A* = (|Aj|, |Aj_i|, |Aj_ 2 |, ■ • •) is the transpose of A, and M is a real matrix of tridiagonal form 



(49) 



M = 



( di c/-i 









\ 



d 3 c 2 
c 2 d 2 c\ 
c\ d\ I 



(50) 



with 



and 



I = J + 1 — m, 



dk — Mfc+m-l 
Cfe = |^/£+m|- 

The largest eigenvalue, xi, of M determines the maximal fidelity through the relation 

2 ' 

To find xi, we set up a recursion relation for the characteristic polynomial of M: 

Qi+i(x) = (d t+ i - x)Qi(x) - c 2 Q t -i(x), 



(51) 



(52) 



(53) 



(54) 
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with the starting values Q-i(x) = and Qa(x) — 1. Eq. |5J resembles the recursion relation of orthogonal polynomials, 
but at first sight the solution does not seems straightforward at all. We thus work out in detail the simplest case for 
which niA = tub = 0. For this particular instance (J54) reads 



Qi+i(x) = —xQi(x) 



-Qi-i(x), 



il 2 - 1 

where we have used the definitions (fff]), (E^) and (E2J). We can rewrite as 



a+i) 



(2/ + l)(2i-l) 
(1 + 1)1 



Qi+i[x) 



= (2l + l)x 



21- 1 
I 



Qi(x) 



- 1[ - Qi-i(x)]. 



(55) 



(56) 



It is now apparent that the terms inside the square brackets can be absorbed into a redefinition of the characteristic 
polynomial through a x-independent change of normalization, namely, 



Qi(x) = (-i) 1 



Pi(x) = (-l) l -^-P l (x). 



(21-1)1! 

This leads us to the recursion relation of the Legendre polynomials: 

(I + l)Pi+i(x) = (21 + l)xPi(x) - lPi-i(x) 



(57) 



(58) 



Working along the same lines, it is easy to convince oneself that the general solution of ( p4[ ) is, up to a normalization 



factor, the Jacobi polynomial P z a,b (x) |] 



where 



= \m B -m A \; b = m B + m A ; 



(59) 



(60) 



and m is defined in (|43|). Note that m can be written simply as m = (a + b)/2. Note also that P l ' is the Legendre 
polynomial P;. 

^,From the result (A12) in Appendix [A] it turns out that the maximal value of the fidelity ( |5^ ) is attained for 
niA — ttib = 0, i.e., precisely the particular case of Legendre polynomials discussed above. Thus, from (|53|) we have 



1 



0,0 
c j+i 



(61) 



where x^ b stands for the largest zero of P„ ,b (x). The fact that itia — tub =0 implies maximal fidelity can be 
translated into physical terms by saying that Alice's states and Bob's projectors must effectively span the largest 
possible Hilbert space. For a fixed choice of itia, not necessarily optimal, the best m B is that for which the Hilbert 
spaces spanned by U(n)\A) and U(n B )\B) coincide, i.e., niA = tib = tn. In this case, the maximal value of the 
fidelity is given by (|53|), with x\ — x j_^_ m , i.e., F = (1 + xJ_^_ m )/2 < F m ax- One reache s the same conclusion if 
m B is fixed and tyia can be adjusted for best results (see discussion in appendix ^ after Eq. A12). 



IV. GENERAL CASE: N SPINS 



We now show that the solution we have obtained in the preceding section is in fact of general validity. Recall that 
in our original problem Alice has N spins. Let us suppose that N is even (N odd will be considered below). As usual, 
Alice constructs her states by rotating a fixed eigenstate of S z . In terms of the irreducible representations of SU(2), 
such state can be written as: 



N/2 I \ AT/2 

\ A )= E \ Y, A ^m A ;a) ; £ £ L4?(m)| 2 = 1. 



(62) 



]=m A a 
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The main difference with the previous example of two equal spins s is that for j < N/2 the irreducible representations 
appear more than once in the Clebsch-Gordan decomposition of (1/2)®". Hence, we label the different occur- 
rences with the index a, which we can view as a new quantum number required to break the degeneracy of Alice's 
system of spins under global rotations. Similarly, the expression for Bob's fixed state \B) is 

N/2 / \ 

\ B ) = E £i#,roB,jS> • (63) 

j=rn B \ J 

However, it is known that equivalent matrix representations 

= (j,m;a\U(fi)\j,m;a) (64) 
are not orthogonal under the group integration, i.e., for a ^ p one has in general 

'dn^$(n)^;{n)^0, (65) 



and the completeness relation J dns O(ns) = I does not hold for the simple choice of projectors 0{tIb) = 
cU(nB)\B)(B\U' ! (ub)- We can circumvent this difficulty by introducing several copies of \B)(B\. A single direc- 
tion (unit vector) Hb is thus associated to 

0(n B ) = U(n B ) [ \B)(B\ + \B')(B'\ + \B")(B"\ + ■ • ■] U\n B ). (66) 

The fixed projectors in the square brackets will be judiciously chosen to eliminate the off-diagonal terms coming from 
the mixing of equivalent representations in the closure relation. The pro jector 0(??b) are explicitly of rank higher 
than one. However, recalling jj), we can view the right hand side of ( |66|) as defining a sum of rank one projectors 
0(ns) + 0'(Hb) + + • • •. The two points of view are equivalent if the averaged fidelity is used as a figure of 

merit. In a suggestive compact notation we can write 

\B){B\ + \B'){B'\ + \B"){B"\+ ••• = \B)-{B\, (67) 

where 

N/2 / \ 

\B)= ]T \ Y^B^j,m B ,f3) , (68) 
j=m B y p J 

and 

B^(B^B'^B"^...). (69) 
Next, we introduce a set of orthonormal vectors {&"}; 

and define the vectors B" as 



fl- = J^±lfty. (71) 



With the above definitions one can easily see that J dns 0(ft B ) = I and, hence, the set of projectors ( pq ) defines a 
POVM. 

The fidelity can be read off from p5|) and it is given by 



N/2 N/2 m-1 

F = 2 + 2T,H Mi(^) 2 + E E *u (^i • ) - \ E E(^) 2 > (^) 



10 



where the phases have been chosen so that Vj, A" and B" are real. In general 6" el 1 , where k must be greater or 
equal than the highest degeneracy of the irreducible representations in the Clebsch-Gordan series of (1/2)®", since 
otherwise ( |70| ) could not be fulfilled. Equation [72] suggests the definition 

J2 A Pl (73) 



which enables us to write 

1 1 



_J2 H \A j \ 2 + £ A^-A^- 1 - J2 l^l 2 - (74) 

j=m j=m+l j=m A 

Using Schwarz inequality we have 

1 N/2 N/2 rn-l 

j=m j =m +l j=m A 

The right hand side is exactly the fidelity (fig) of the preceding section with the substitution 



^^^W = W (76) 

This equation shows that the existence of several equivalent representations in the Clebsch-Gordan decomposition 
of Alice's Hilbert space cannot be used to increase the value of the fidelity already obtained in section III. The 
equality holds when all vectors Aj are parallel, in which case we recover J45|). The square root on the right hand 
side of (|7q) plays the role of an effective component of \A) on the Hilbert space of a single irreducible representation 
j. The specific ways \A) projects on each one of the equivalent representations are of no relevance, provided Aj 
do not change. As far as the fidelity is concerned, all them are equivalent to taking a state \A) that belongs to 
N/2 (AT/2 — 1) © (N/2 — 2) • • • (no duplications), with the corresponding components given by Aj. 

As we have just seen, the maximal fidelity can be achieved from a code state containing only one of each irreducible 
representations. This type of states are formally the same as those considered in the simplified example of two equal 



spins si — s 2 — s studied in section III. for which s (g) s — J © (J — 1) © • • • © 0, with J = 2s = N/2. The problem 



of an even number of spins is thus completely solved: according to 153) the maximal fidelity is given by 

1 + 

F N = ^±1, for N even, (77) 

where 2^/ 2 +i 1S the largest zero of the (Legendre) polynomial Pn/2+i( x ) = ^72+1 0*0 • 



For an odd number of spins we can proceed as in section [II but considering now states with two different spins: 
si = s, S2 = s + 1/2. The corresponding Clebsch-Gordan decomposition is also non-degenerate: s g) (s + 1/2) = 
J © (J - 1) © ••• © 1/2, with J = 2s + 1/2 = N/2. The results from @ to @ are still valid (for the value of 
J we have just specified). The optimal values of and mg are again the minimal ones: ttla — tub = 1/2. The 
maximal fidelity is 

1 + x ' 1 

F N = X JLp±l/l i for N odd; ( 78 ) 

where 2^/2+1/2 s t an ds for the largest zero of the Jacobi polynomial P^j 2 +i/2^ x ^ This completes the solution of the 
general problem. 

It is physically obvious that the larger the number of spins Alice can use the better she should be able to encode 
71. One thus expects that the maximal fidelity should increase monotonously with N. It is interesting to obtain this 
result from the properties of the zeroes of the Jacobi polynomials. For an even number of spins, N = 2n — 2, the 
corresponding zero is x%°, whereas for TV + 1 it is a;^' 1 , and x^_ x for N — 1. Proving that Fjy_i < Fjq < Fjv+i 
amounts to showing that 



< x n ' < a; n ' , (79) 
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N 


1 


2 


3 


4 


5 


6 


7 


F N 


2 
3 


3+v^ 
6 


6+Ve 

10 


5+v^5 
10 


.9114 


.9306 


.9429 



TABLE I: Maximal fidelities as a function of the number of spins. 



but this is just a particular case of (A9) for a = and 6 = 1. 

Not only the optimal strategy Alice can devise with N spins leads to a fidelity larger than .F/v-i- She can also use 
non-optimal ones and still exceed .Fjv-i- E.g., for AT = 4, the choice tua = m-B = 1, which is non-optimal, gives a 
fidelity F = (10 + yl0)/15 > (6 + v6)/10 = i^s. This is also a trivial consequence of (A9) as in this case one has 



0.2 



> x. 



0,1 



In physical terms, this is telling us that the dimension of the Hilbert space spanned by U(n)\A) and 
U(nB)\B) when N = 4 and tua = TOb = 1 (including equivalent spin representations only once) is still larger than 
the maximal available dimension for AT = 3. 



V. DISCUSSION AND OUTLOOK 



In this paper we have addressed the problem of optimizing strategies for encoding and decoding directions on the 
quantum states of a system of N spins. We have restricted ourselves to states that point along a definite direction in 
an intrinsic way, namely, to eigenstates of ft ■ S. This case is of great interest since no prior knowledge of any sender's 
(Alice's) reference state or frame by the recipient (Bob) is required at all for a viable transfer of the information. We 
have optimized both Alice's states and Bob's measurements. Our results are summarized in ( fT7| ) and (|78|), where we 
give the maximal averaged fidelities Fjy. Interestingly enough, these results can be written in terms of the largest 
zeroes of the Jacobi polynomial, which are known to play an important role in angular momentum theory and are 
intimately related to the matrix representations of SU(2). The states that lead to the maximal fidelities are among 
those that have the smallest (non- negative) value of n • S, namely, to = for N even and to = 1/2 for N odd, but 
still span the largest Hilbert space under rotations. 

We display the values of the maximal fidelity for A~ up to seven in table Sfor illustrational purposes. It shows, 
e.g., that the optimal encoding with three spins (m = 1/2) gives F3 = (6 + V6)/10 ~ 0.845, which is already larger 
than the corresponding maximal value for four parallel spins (to = 2): F — 5/6 ~ .833 0. This illustrates a general 
feature: the optimal strate gies d iscussed here lead to fidelities that increase with A" much faster than that of sending 



parallel spins. In fact, Eq. A13 shows that Fjy approaches unity quadratically in the number of spins, namely 



F N ~l-£-, (80) 

where £ ~ 2.4 is the first zero of the Bessel function Jo(x). In contrast, if parallel spins are used the maximal fidelity 
approaches unity only linearly, F ~ 1 — 1/N. 

This can be understood in terms of the dimension d of the Hilbert space used effectively in each case, which is a 
direct sum of the Hilbert spaces of the irreducible representations of SU(2) involved. Here effectively means 'non- 
rcdundantly', thus equivalent representations count only once. Encoding with A^ parallel spins uses only the Hilbert 
space of the representation N/2, whose dimension is d = N + 1, whereas our optimal strategy uses a much larger 
Hilbert space, with d = (N/2 + l) 2 for N even and d = (N/2 + l) 2 - 1/4 for N odd; in both cases d - A^ 2 . We are 
led to the conclusion that the fidelity as a function of d tends to unity as 

F~l-- d , (81) 

where a is of order one and depends on the particular strategy. 

Improvements on the approach discussed in this paper can only come from encoding and decoding procedures that 
make extensive use of the available Hilbert space, namely, strategies that use the redundant equivalent representations. 
In fj] we presented a strategy for which the maximal fidelity approaches unity exponentially in the number of spins, 
i.e., F ~ 1 — 2~ N . We argue there that this encoding is likely to lead to the maximal fidelity one can possibly achieve 
with N spins, since it makes effective use of the whole Hilbert space of the system, for which d = 2 N (thus, Eq. |8l] 
also holds in this case). The corresponding encoding process, however, involves complicated unitary operations and, 
moreover, it seems to require that Alice and Bob share a common reference frame [ fL3| . 

We have obtained our general results using continuous POVMs, but finite ones can also be designed. For A^ parallel 
spins (m,A — rns — N/2), a general recipe for finite optimal POVMs exists H, and minimal versions for up to 
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N = 7 can be found in Q. The unit vectors ri r associated to the outcomes of these POVMs are the vertices of 
certain polyhcdra inscribed in the unit sphere. For TV < 7 we have explicitly verified that these very same polyhedra 
can be use to design finite optimal POVMs for any value of ttia — rns < N/2. Moreover, the minimal POVMs 
of (J] remain minimal for the states considered here. We have discussed this issue in detail for N = 2 in section II. 
For TV = 3 the polyhedron corresponding to the minimal POVM is the octahedron M. One can easily verify that 
O r = U(n r )\B}(B\W (n r ), where \B) is given in ( |39| ) with tub = 1/2, 3/2, fulfill the completeness condition (||). We 
hence believe that the discretization of a continuous POVM is a geometrical problem, i.e., it seems to be independent 
of the states \B). 

The optimal states, \A), can be easily computed from the matrix M in (|5C|), as they are the eigenvectors corre- 
sponding to the maximal eigenvalue. Recall that for N — 2 one obtains the one-parameter family of states (^J) which 
includes the product states | ||), | []). For N > 2, product states of the type | Tittl ' ' ') do not seem to be optimal. 
Consider, e.g., N = 4. The optimal eigenvector of M is 




| A) = ^-|2, 0) + e *» -j= 1 1, 0) + e**> ^ - 10, 0) , (82) 

which is clearly not a product state of the individual spins for any choice of the phases |4j . One could argue that this 
solution is not entirely general because the Clebsch-Gordan series of (1/2)® 4 contains the representation 1 three times 
and twice, whereas in ( p2[ ) they appear only once. However, any optimal state has the same 'effective' components 
Aj (see eqs. fn||76|), which can be read off from (|82|): 

i2 = TP il = 7r (83> 

Note now that any product state with m — (two spins up and two spins down), e.g., | ttll); I TIIT)i nas an 'effective 



Clebsch-Gordan decomposition given by A2 = A\ = Aq = l/\/3, which do not coincide with (83). Therefore, these 
product states cannot be optimal. Nevertheless, they lead to a maximal fidelity F — (15 + 5\/2 + 2v / 5)/30 w 0.885, 
which is remarkably close to F4 w 0.887. This is likely to be the case for arbitrary N. These issues are currently 
under investigation. 
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APPENDIX A: 

In this appendix we collect the mathematical properties of the Jacobi polynomials P^ b (x) that we use in the text. 
We are concerned only with integer values of a and b such that b > a > 0. Further properties can be found in p5Q 
and |i . 

For fixed a and b, {P£' b (x)} is a set of orthogonal polynomials, where n labels the degree of each polynomial in the 
set. A convenient definition can be stated in terms of their Rodrigues formula: 

P%\x) = - *r a (l + *r b ^T [(1 - z) n+a (l + xT +b ] ■ (Al) 

From (Al) it follows the recursion relation: 

xP»> b {x) = arjtf^x) + (3 n PZ< b (x) + 7„P:i\(z), (A2) 



with 



A. 



2(n + l)(n + a + b+ 1) 
(2n + a + b + l)(2n + a + b + 2) ' 

b 2 ~a 2 

(2n + a + b)(2n + a + b + 2) ' 



2(n + a)(n + b) 

ln (2n + a + b)(2n + a + b+l)' 1 ' 
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It also follows the differentiation formula 

d P^ b (x) _n + a + b+l DO+ i,6+i 



dx 



-P 



(x). 



(A4) 



The normalization is chosen so that the coefficient A n of the highest power of P^ b {x) = A n x n + B n x" 



is 



A r . 



F(2n 



2™n!r(n- 



The following two relations can also be obtained from the definition Al 



(2n + a + b)P 1 a /- 1 (x) = {n + a + b)P r a l ' b (x) + (n + a)P^ 1 (x) 



a,b 



(n + b + a + 1) 



1 



-P 



a,b+l 



(x) = (n + l)P^-\x)+bP^(x) 



ya,b/ 



(A5) 

(A6) 
(A7) 



Let us recall some basic facts about the zeros of orthogonal polynomials: i) any orthogonal polynomial, P n , of order 
n has n real simple zeros. For Jacobi polynomials these zeros lie in the interval (—1, 1); ii) the zeros of P n and P n +\ 
are interlaced; Hi) for x greater th an t he largest zero, the polynomial is a monotonously increasing function (if the 
polynomial is normalized as in Eq. A5, where A n > 0). In particular, P n (x) must be positive in this region . 

Now we can prove the results needed in the text. As in there, we denote by x^ b the largest zero of the polynomial 
P%' b (x). Let us start by showing that 



a+1,6+1 a ,b 



(A8) 



From property Hi above it follows that the left hand side of (A4) is manifestly positive for x > x^ b . Hence, so it is 
the right hand side. We conclude that x^t^ +1 cannot belong to this region and (AS) follows ■ 
Next, we prove the inequality 

(A9) 



We evaluate (A6) at 



and use properties ii (=>- x n-i < x n b ) and Hi, which imply that Pn-i( x n' b ) > 0, to 



a , b 



obtain that P r a l ' b ~ 1 (xf l - b ) > 0. We repeat the process for x = ar^ x and conclude that P%' b ~ 1 ( x % b i) < 0. Hence 
pa,fc-i nag a zer0 j n m ^ erva i ( x n-ii x n' b )- This is necessarily the largest zero a^ ,fc_1 since, according to (A6) and 
properties ii and Hi, P^' b ~ 1 (x) > for x > x^ b . Thus ( |A9| ) follows ■ 
The inequality 



„a,b+l ^ _o,6-l 
x n ^ x n+1 



(A10) 



can be proven as follows. Evaluate (A7) at x — x^ b+1 so that the left hand side of this equation is zero. The second 
inequality in (A9) and property Hi imply that P®' 6 (x£' + ) > 0. Hence the first term on the right hand side of (A7) 
must be negative, i.e., Pn+i 1 ( x n' b+1 ) < 0, and ( AlOj ) follows immediately, since otherwise property Hi would not hold 

for p^-'m 



For two given integers I, m consider now the following set of zeros 



CI 



r m 

i x l- 



m —m ,m +m 



777 < Tfl < TTl < 



(All) 



We want to prove that 



max CL 



"i-m ■ 



(A12) 



According to (lAa), decre asing m" by one le ads u s to a larger zero. The maximum is then in the subset {x 



m' < I}. The inequality QA1Q ) now implies (A12) ■ 

Finally, we give the large n (asymptotic) behavior of x°^ b jE 



0,2m' 
l—m' 



m < 



£ 2 l 

2n z nr 



(A13) 
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where £ a is the first zero of the Bessel function J a (x). For a = 0, which is relevant for our discussion in section [V], we 
also give the subleading term: 

4< b = l-^(l-— )+0(^), (A14) 

where 

& = i = 2.405. (A15) 
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